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Abstract: By employing the refinable matrix of the scaling function, the method for constructing 
a pair of tight framelets with differential relations on the interval is given in this paper. 
In which the tight framelets is defined in the sense of Chui’s definition. Moreover, an 
explicit example is constructed by using B-spline. This result can be applied to solve some 
differential equations. 
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1 Introduction 


Divergence free wavelets have some applications for the analysis and numerical simulation 
of incompressible flows typically modeled by the Stokes system or by the incompressible Navier- 
Stokes equations!!]. In order to find the kind of wavelets over bounded domain, two families 
of wavelets with differential relations on the interval are needed firstly. Although there exist 
many wavelets on the interval?~4], none of them has differential relations. Lakey and Pereyral®! 
constructed wavelets on the interval with differential relations by using the multiwavelets of 
L2(R), which was introduced by Hardin and Marasovich!§l. More precisely, by truncating the 
scaling functions and the wavelets of L2(R), they obtained a pair of biorthogonal MRA’s and 
wavelets for L2[0,1]. Then they found another pair of MRA’s and wavelets by smoothing 
and roughing method. It turns out that the two pairs of wavelets and scaling functions have 
differential relations. However, the corresponding divergence free multiwavelets have too many 
generators and duals, which cause problems in applications. 

Framelet is a widely accepted tool in signal and image processing as well as in numerical 
simulation. In particular, they perform better effects than wavelets in many fields of image 
processing. Ron and Shen!) introduced framelets MRA of L2(R4%) with the setting of shift- 
invariant subspaces, and give the filter bank algorithms of the pyramidal decomposition and 
reconstruction. Since the construction of tight framelets over bounded domain is quite different 
from that over unbounded domain, Chui et al! developed a general theory for construction of 
tight framelets over a bounded interval using univariate splines. Due to the theory, this paper 


constructs framelets with differential relations on the interval. 
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2 Tight Framelets on the interval 


In this section, we introduce some notations and preliminary results, which will be used 
later on. To construct tight frames, we need a sequence of nested subspaces {Vj}jen in Lo(J) 


which satisfies 


Vico CVj C+ => L(I), and (J Vj = (1), 
j21 
where V; is linearly spanned by ¢;,1,°-- ,%j,m;- Although this sequence does not have both 
translation and dilation invariant properties, it generates a multiresolution analysis (MRA) over 
bounded domain I. Let ©; denote the column vector [6j,1, + ims and Mj = {1, , m5}. 
Then there exists a matrix Pj of size mj x mj41 such that ©; = Pj®j+ı due to V; C Vj+1. 
The matrix P; is called a refinable matrix. A function family {®;};>1 is said to be locally 
supported, if the sequence 
h(®;) := mor length(supp¢;,x) 


converges to zero. 

The next two definitions!®! are also needed in our paper. 

Definition 2.1 Let ®; be a finite family with cardinality m; in Lo(I ). For any symmetric 
positive semi-definite (spsd) matrix S; = [sË )] k lem; the quadratic form T} is defined by 


T= Worren a ens Je LO: (1) 


and the corresponding kernel 


Ks, (z,y) := > sy) j,k (2) 5,€(Y)- 
k,lem; 
Definition 2.2 Assume that {®;};>1 is a locally supported family and Sı is an spsd 
matrix, which defines the quadratic form T, in (1). Then the family {Vj}j>1 = {Qj®j4i}j>1 
constitutes an MRA tight framelets of L2(Z) with respect to $1, if 


Tif +> Yo Mf dja)? = A, VF © Lol). (2) 
JRL KEN; 

We find that the ground level S is relevant to the order of vanishing moments of the 
framelets Y; g. The following theorem tells us when V;,;, constitutes tight framelets. 

Theorem 2.1!! Let {®;}j>1 be a locally supported family and Sı an spsd matrix such 
that |T; fl] < || f||? for all f € Lo(Z). Then {Vj}j;>1 = {Q;® 541} forms an MRA tight framelets 
with respect to S,, if and only if there exist spsd matrices S; of dimensions m; x mj, j 2 1, 
such that the following conditions hold: 

(i) The quadratic forms T; in (1) satisfy 


jim Tf = (IFIP, f eLa; (3) 
(ii) For each j > 1, Qj, S; and $;41 satisfy the identity 


S541 = P? S;Pj = QF Q;. (4) 


NO. 2 Zhuang Zhitao, Luo Dang: Tight Framelets with Differential Relations on the Interval 281 





Remark 2.1 The condition 


lim Ks, (x, y)ldy = 0 (5) 
|z—y|>e 
for any € > 0 is sufficient for property (i) in Theorem 2.1. If the matrices S; have a fixed 
maximal bandwidth r > 0 and {®;}j51 is locally supported, then (5) holds, since the integral 
in (5) is zero for sufficient large j. 

Remark 2.2 If S; = I; for j € N, then (4) reduces to Im,,, — P7 P; = Q7Q;. In this 
situation, Lai and Naml?) gave a simple construction of tight framelets on the interval. However, 
their framelets don’t have any vanishing moments. 


3 Differential relation 


This part is devoted to construct two families of tight framelets with differential rela- 
tions on the interval. Let {V jz and {V; }j>ı be two nested subspaces of La(/), the 
corresponding generators are oF and ®;, respectively. We are interested in the differen- 


tial relation ie = M;®;, where M; is an mf x my matrix, ®f = [b)1,--+ ,@jm;]” and 
’ / : 
op = OR f y Olm | Since OF = PHS, ‘hat 
+ tpt” 
Dt = PTO, | = P Mj 85,1. 


_ ie Ab aban aed , = cree z 
On the other hand, ®; = P; ®;,; implies oF = Mj; = MjP; ®;7,,. Suppose {54,4 }kemy 4: 


are linearly independent, then we get 
Mj Py = P} Mjr. (6) 


Under this assumption, the Gramian matrix Tj = [(05,, 0; ¢)lk, tem; is symmetric positive 
definite, and its dual basis S7 is given by the function vector S- j = (6 by alkem; = T3 are It 
is well known that 
2 os That = 
MAI = [Hore (Foe 
Now, we give a simple method to construct two families of tight framelets with differential 
relations on the interval. 

Theorem 3.1 ee A 1 and {®; }j>1 be two locally supported families with differen- 
tial relation oy = 7 {UP = Q7 ® tak be an MRA tight framelets with respect to an 
spsd matrix ST. ee 

= ae 1 5 
W =Q bjn = 523 Min P54 
for j > 1, where b is a positive number in R such that 
-1 Il r 
=M 
is an spsd matrix. If the matrices 57 := aM; TST Mj; have fixed maximal bandwidth r > 0, 
then {V; }j;>1 constitutes an MRA tight framelets of L2(I) with respect to Sy and př = = bW;. 
Proof Since S} is an spsd matrix, S; has the same property. Note that 


KF, CRSA A K, $1) | keM, 


ry Si M, (7) 
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is the norm of the orthogonal projection of f onto Vy . Then for f € L(I), the condition (7) 
ensures that 


if = KF, daaa [F binl kem, S + (Cf, $i, Pieris KS, 1, pen = (FU? 


and T; satisfies (i) of Theorem 2.1 due to the Remark 2.1. Now, we only need to prove S7 
satisfies (ii) of Theorem 2.1. In fact, 


seilh 
ma P SP = F St) Mya1 — 


aM. =z T TMI S} M;P7 


b 


T 
j L (MZ 5} a Mjs1 — MI PE St Pt M341) 


1 T 

pein (Sp PY SP PP) Mian 

1 T ad yen 

b2 Ming} Ql Mjn = Q; Q; ’ 

where the second equality holds by using the fact (6) and the forth equality due to (ii) of 
Theorem 2.1. Thus {®; }j>1 constitutes an MRA tight framelets of L2(J) with respect to S7. 


Moreover, we have 
1 + a, ae as = 
vy = = =O) = = Qj Mjn = bQ; ip = bY. 


This completes the proof. 

For the efficiency and simplicity of computation, B-spline is usually used for constructing 
wavelet functions. Next, we apply B-spline to Theorem 3.1 in order to find the desired tight 
framelets on the interval. 


Let us recall the scaling relation of B-spline ¢” for m > 2 (see [10]) 


p(s) = X g'o (2x — k), 


kez 


where 
27MH (7), for O0<k <m; 


Ck = 
0, otherwise. 


Consider B-spline function ¢” of order m whose dyadic translations are restricted into interval 
(0, m], i.e., 6 (29 - —k)|fomj. Assume 
FC) = 2i-1 em (ai? . —k)|fomj, and Vj" = {97k : 1< k< mi}, 


then the family of nested sequence {V;” : j € Z4} forms an MRA of L2([0, m]). Let 


Pr := [P70 ie 


J jl? ETM; 


Then their exists a mj x mj41 matrix P™ such that 7 = Prd , for each j € N. In [9], Lai 
and Nam gave a method to construct tight framelets in this situation. We show the quadratic 


B-spline and its tight framelets of the ground level in Figure 1. 
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; , = _ ? spay 
Since (™) (z) = =! (x) — -Hx — 1), (P7) := M; with 
—1 
1 -1l 
; 1 -l 
M; = 27 
—1 
1 
m 5 x(m—1); 
Define SP- := gs MF, ,Mj41, where b € R+ is chosen such that 'y'~* — ST is an spsd matrix. 


Obviously, S has a fixed maximal bandwidth 2 due to the form of M;. Hence, we can find 


the tight framelets 
= 1 = 
vy = gY Mmda 


By Theorem 3.1, y := b47. Take m = 3 and b = 2, then 


Boe. ā 4 2 
15 30 15 15 

saa ia SS i 

p27} g2% = 30.15 30 «15 
: g 4 1343 B B 
15 30 15 30 

oe 4 B 3 

15 15 30 15 


is an spsd matrix. Hence, y3 and v5 are tight framelets with differential relation p? 
We show linear B-spline and their tight framelets of the ground level in Figure 2. 


Figure 1: B-spline ©} (left) and tight framelets YÌ? (right) 





— 9y3- 
=20};. 
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Figure 2: B-spline ©? (left) and tight framelets YY“ (right) 
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